################################################################################################
# Check for symmetric meaning of "Ban" and "Permit"
options(warn=-1)
#data$rel_posban <- droplevels(data$rel_posban)
#data$rel_posper <- droplevels(data$rel_posper)


cbind(round(prop.table(table(data$rel_posban)), 2))
cbind(round(rev(prop.table(table(data$rel_posper))), 2))


ban <- as.vector(prop.table(table(data$rel_posban)))
perm <- as.vector(rev(prop.table(table(data$rel_posper))))
p.chi <- chisq.test(ban, perm)

ban <- round(ban, 2)
perm <- round(perm, 2)

#ban <- as.vector(table(data$rel_posban))
#perm <- as.vector(rev(table(data$rel_posper)))
#chisq.test(ban, perm)
#options(warn=0)

